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ABSTRACT 

We study structure formation in a phenomenological model of modified gravity which in- 
terpolates between ACDM and phenomenological DGP-gravity. Generalisation of spherical 
collapse by using the Birkhoff-theorem along with the modified growth equation shows that 
the overdensity for spherical collapse S c in these models is significantly lowered compared 
to ACDM, leading to enhanced number densities of massive clusters and enhanced cluster 
merging rates. We find that 6 c (z) is well fitted by a function of the form S c (z) = a — b exp(-cz). 
We examine the sensitivity of PLANCK's and SPT's Sunyaev-Zel'dovich survey to constrain 
the modified gravity parameterisation and find that these experiments can easily distinguish 
between models with a cosmological constant and modified gravity, if prior constraints from 
CMB temperature and polarisation anisotropics are included. 
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The late-time accelerated expansion of the universe, observed in 
the magnitude redshift relation of type la supernovae, can be ex- 
plained by three very different theories, either by assuming that 
the expansion is due to the cosmological constant A, or due to 
dark energy as an additional fluid, or due to a modification in the 
gravitational law. Although there has been rapid progress in devel- 
oping modified theories of gravity to explain the late-time accel- 
erated expansion of the universe, none of the models is yet con- 
vincin g. The Dvali-Gabadadze-Porrati (DGP) model jDvali et all 
2000) for example, which is based on a five-dimensional brane 
world set-up, offers a simple explanation of the accelerated ex- 
pansion of the universe from 5d physics without dark energy, but 
it has been show n that the model suffers from a ghost instability 
iKovamal J2007ah . Simple modifications of the 4d Einstein-Hilbert 
action such as /(if) gravit y violate s olar system constraints un- 
less f(R ) is fine-tuned dChibal |2003|). Structure fo rmation with 
DGP-gravity was investigated by Lue et al.l I l2004al |b|) , and con- 
straint s on DGP gr avity have been deriv ed (CMB: ISong et all 
d2007l). supernovae: | Rydbeck et al] d2007l). large-sc ale structure: 
iMaartens & Maierottol d2006t): iKovama & Maartenj {2006); a re- 
view is provided bv lKovamafeoOVbh ). 

However, we have learnt general features of modified grav- 
ity from these attempts: Firstly, any modification of general rela- 
tivity (GR) on cosmological (horizon) scales results in modifica- 
tions of GR on sub-horizon scales as well, which influences the 
growth of perturbations. Whereas in cosmologies with homoge- 
neous dark energy the growth function is uniquely determined by 
the Hubble-function, modified gravity models generically alter the 
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growth function and it becomes impossible to reconcile measure- 
ments of the expansion and growth rates resulting from a modi- 
fied gravity model with a dark energy model with standard grav- 
ity. This leads to the interesting possibility of distinguishing the 
modified gravity models from dark energy models using probes of 
the cosmological large-scale structure such as the CMB or the inte- 
grated Sachs- Wolfe effect. Secondly, any successful modified grav- 
ity models should satisfy the stringent constraints coming from the 
local tests of gravity. This effectively means that the theory must ap- 
proach GR on small scales. The DGP-model naturally recovers GR 
on small scales without introducing any additional phenomena due 
to the strong coupling. An example of this is /(/f)-gravity, where a 
very contrived form of the f(R )-function (the so-called Cham eleon 
mechanism) is necessary to recover GR dHu & SawickHl2007t) . 

An interesting probe of the growth of structure is provided by 
weak lensing, but the non-linear transition to GR is important for 
the correct calculation of the weak lensing signal. The common ap- 
proach in the literature is to parameterise the linear growth rate and 
to use the non-linear mapping of the linear matter power spectrum 
developed in GR to calculate nonlinear growth and the nonlinear 
lensing signal. This is effectively equivalent to assuming that grav- 
ity is modified even on small scales which would contradict the 
local gravity experiments. Recently, attempts ha ve been made to 
carry out w-body simulation s in modified gravity (Sta benau & Jainl 
2006; Laszlo & Bea"nll2007t) . but the mechanism to recover GR due 
to non-linear dynamics was not properly implemented. An interest- 
ing phenomenological attempt was ma de to address this non -linear 
transition to GR using the halo-model dHu & Saw icki 2007), but it 
is still an open question to calculate the non-linear power spectrum. 

In this paper we investigate constraints on modified gravity 
as an alternative to dark energy dDvali & Turnerll2003l : iKovamj 
2006) from structure formation and from the abundance of clus- 
ters of galaxies by solving the homogeneous growth and the spher- 
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ical collapse equations. As the model for modified gravity we em- 
ploy a phenomenological model which is able to interpolate be- 
tween ACDM and DGP-gravity, and predict the number density 
of clusters of galaxies in a Press-Schechter approach. Due to the 
spherical symmetry, it is possible to address the non-linear dynam- 
ics of t he gravity in a trac t able way. Similarly to the case of dark 
energy I lLahav et alJll99ll : iMultamaki et~afll2003l) . cluster counts 
are a promising tool for investigating modified gravity as the col- 
lapse threshold for objects forming at late times is lowered, which 
gives rise to an enhanced spatial density of clusters. The influ- 
ence of dark energy on cluster formation was extensively studied 
dWang & Steinhardllll998l : |Percivalll2005l ; iBartelmann et alj|2006t 
iMaor & Lahavll2005t to name but a few), and we aim to forecast 
constraints on modified gravity from cluster count experiments. 

This paper is structured as follows: The theory is developed 
in Sect. [2] where our model for modification of gravity is outlined. 
The results are compiled in Sect. [3] for homogeneous cosmology, 
in Sect. [4] concerning the solution of the spherical collapse equa- 
tions, and in Sect.[5]fhe impact of modified gravity on cluster num- 
ber counts and merging rates is discussed. Constraints on cluster 
formation with modified gravity resulting from future Sunyaev- 
Zel'dovich surveys are given in Sect. [6] A summary in Sect.|7]con- 
cludes the paper. Parameter choices for the fiducial cosmological 
model are based on the spatially flat ACDM cosmology with adi- 
abatic initial perturbations: Specific choices for the cosmological 
parameters are Sl„, = 0.25, Ho = 100/ikm s~' Mpc~' with h = 0.72, 
a„ = 0.04, n s = 1 and o" 8 = 0.8. 



2 STRUCTURE FORMATION IN MODIFIED GRAVITY 
2.1 Phenomenological DGP-gravity 

A lack of satisfying theoretical models for modified gravity en- 
forces us to make phenomenological approaches. We assume that 
the modified Friedmann equation (normalised to H ) is given by 



H 2 (a) = {l-Q. m )H a (a) + 



"IB 



(1) 



where a smoothl y interpolates between ACDM (a = 0) and DGP- 
gravity (a = 1) dDvali & T urner 2003b. We restrict ourselves to 
linear structure formation as well as on object formation by spher- 
ical collapse, where we incorporate modified gravity in the growth 
equation and use the Birkhoff-fheorem for describing the spherical 
collapse dynamics instead of solving the structure formation equa- 
tions TorJML DGP-gravity. 

lKovamal (2006) show that a covariant effective theory that re- 
produces eqn. (TJ can be constructed for a = 1/n where n is an 
integer number. Although the Bianchi-identity would constrain the 
structure of the theory, we need to supply additional information 
in order to derive the equations that govern the evolution of struc- 
ture. In this paper, we choose the Birkhoff's theorem as a guidance 
to study the evolution of structure which ascertains that the metric 
observed by a test particle outside a spherically symmetric matter 
distribution is equivalent to that of a point source of t he same mass 
located at the centre of the sphere dLue et ai]|2004bl) . Recently, it 
has been argued that the violation of Birkhoff' s theorem poses se- 
vere conceptual and computational difficulties dDai et al1l2007l) : If 
the gravitational force in the interior of a spherically symmetric 
matter distribution is non-zero, one would need to include the de- 
tails of the host structure a galaxy is residing in when determin- 
ing gravitational forces inside a galaxy. We should note that the 
DGP-gravity does not respect Birkhoff's theorem. Nevertheless, a 



Birkhoff-like cancellation of the gravitational acceleration inside a 
spherical mass is enforced from the five-dimensional physics. Thus, 
the Birkhoff's theorem is surely not the unique way to avoid the 
above described conceptual difficulties, but it should give reason- 
able approximation to the behaviour of gravity in viable modified 
gravity models. In fact it was shown that the difference between us- 
ing Birkhoff's theorem and the DGP-like gravity is small in general 
cosmolog ies described b y eqn. dTJ concerning linear growth rate of 
structure dKovamal l 2006h . 



2.2 Birkhoff's theorem 

In this section, we explain how the Birkhoff's theorem can be used 
to describe the growth of overdensities. The modified Friedmann 
equation can be rewritten in the form 

SjtG P 
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(2) 



Taking the time derivative of the Friedmann equation, we obtain 

3 



g(fo)-^og'©>) 



(3) 



where g'(£) = dg(^)/d^. Let us consider a spherically symmetric 
collapsing dust shell with radius r(f). If Birkhoff's theorem holds, 
the dynamics of r(t) can be determined as 
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The overdensity is defined by 

r Pshell - P 



(4) 



(5) 



(6) 



Initially, the shell is expanding due to the expansion of the Universe 
such that initially r(t) = a(t)ro and p s heii = p. The conservation of 
mass in the shell means that r(t) and d are related by the condition 
r 3 (0Psheii = a3r oP and it is possible to derive the evolution for 6 
from the equation for r(t). The linearised evolution equation for 
S <K 1 is then given by 



5 + 2HS = 4^G[g'(fo) + 3f g"(fr)]p& 



(7) 



Eqns. 10 and are the basic equations that describe the growth 
of structure in this model. 



3 HOMOGENEOUS COSMOLOGY 

The Hubble function H{a) follows from the Friedmann equation 
numerically by solving eqn. |QJ as an algebraic equation with the 
parameter a. All higher derivatives follow analytically by differ- 
entiating both sides of the Friedmann equation and isolating the 
respective derivative, e.g. for dH/da one obtains: 



dH 3£X„ / ^1 \- 

-r- = -f{a(l-n m )H-'-2H(a)) 



(8) 



da a 

For visualisation of the Hubble function and its derivatives, the re 
spective SCDM scaling is divided out: 

d"H 
da" 



H n (d) : 



(9) 
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Figure 1. Hubble function H(a) and its first and second derivatives, in the 

2n+3 

notation H„ = a 2 d"H/da", n = 0, 1,2, for a = (corresponding to 
ACDM, dash-dotted), a = 1 (DGP-gravity, solid) and intermediate values 
in steps of Aa = 0.2 (dashed). 
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scale factor a 

Figure 2. Evolution dV/dz of comoving volume V with redshift z, as a func- 
tion of scale factor a, for a = (corresponding to ACDM, dash-dotted), 
a = 1 (DGP-gravity, solid) and intermediate values in steps of Aa- = 0.2 
(dashed). 



The results for H„(a) are depicted in Fig. [Tj The effect of modi- 
fying gravity according to eqn. Q} is a more gradual increase of 
the Hubble function Ho(a) compared to ACDM, which is reflected 
by smaller values of the derivative i?i(a) and a smaller curvature 
H 2 (a). 

With the Hubble function H(a), the angular diameter distance 
in units of the Hubble distance c/Hq can be defined: 



d A (a) 



•J a 



da ■ 



1 



a 2 H(a) 



(10) 



which yields the evolution of comoving volume with redshift, in 
units of the Hubble distance cubed (c/H ) 3 , for 4n solid angle: 



dV _ tP A (a) 
— (a) = 4k— . 

dz a 2 H(a) 



(11) 



The evolution dV/dz of comoving volume as a function of scale 
factor a is given in Fig. [2] The models vary most at a - 0.3, where 
they exhibit a variation of almost 20%. The ACDM model shows 
the most dramatic increase of comoving volume as a function of 
redshift, which is due to the rather fast transition between the matter 
dominated and A-dominated phases. 

The evolution of cosmic time with redshift, in units of the 
Hubble time l/H , follows directly from the definition of the Hub- 
ble function as the logarithmic derivative of scale factor a with cos- 
mic time t, H = dlna/df, and the relation between redshift z and 
scale factor a, a = 1/(1 + z): 

dt a 

-Ha) = 77TT ( 12 > 
dz H(a) 

The evolution of cosmic time df/dz as a function of scale factor a 
is plotted in Fig. [3] The influence of modified gravity seems to be 
minor, amounting to a few percent at a = 0.6. Similar to the case of 
dK/dz, the modified gravity models show a slower evolution df/dz 
compared to ACDM. Appendix [A] contains a detailed comparison 
of modified gravity models with dark energy models with identical 
expansion histories, focusing on the effective dark energy equation 
of state w(a). 

The last point is the growth equation, which follows from the 
linearised structure formation equations and describes the linear 
and homogeneous growth of the cosmic density field with time, 




scale factor a 

Figure 3. Evolution df/dz of cosmic time f with redshift z, as a function of 
scale factor a, for a = (corresponding to ACDM, dash-dotted), a = 1 
(DGP-gravity, solid) and intermediate values in steps of Aa- = 0.2 (dashed). 



S(x,a) = D + (a)6(x,a = 1). The growth equation with the scale 
factor a as the time variable reads (c.f. eqn[7](: 



d 2 1 

7T D + + - 
da z a 



3 + 



din H(a) 



din a 



^-D + = —G(d)a m (a)D + (a). (13) 
da Za z 



The modificat ion in the gravita tional law is in corporated in the 
function G(a) dLue et alJ2004ah : 



G(a) 



2a 4 
30T 



dH dH 

3H(a) + a\ 

da \ da 



+ aH(a) 



d 2 H 

da 2 



(14) 



and is reflected in the epoch-dependence of the matter density pa- 
rameter f2,„(a), 



Q. m (a) 



(15) 



s H 2 (a)' 

Fig. [4] shows the growth function D + (a) resulting from numerical 
integration with initial conditions D + = and d£> + /da = 1 at 
a = 0, with a as a free parameter. The normalisation is chosen 
such that D + = 1 today. Common to modified gravity models is 
a faster growth of structure at early times, followed by a slower 
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scale factor a 

Figure 4. Growth function D + (a), for a = (corresponding to ACDM, 
dash-dotted), a = 1 (DGP-modified gravity, solid) and intermediate values 
in steps of Aa = 0.2 (dashed). Additionally, the growth function D + (a) = a 
for the SCDM cosmology was plotted. 



structure growth from the point on where the modification of grav- 
ity becomes important. The largest spread between the models with 
different a can be observed in the vicinity of a = 0.4, where the 
growth functions differ by almost 10%. Appendix IB"1 provides a de- 
tailed discussion about the magnitude of the terms governing the 
evolution of the growth equation eqn. J 1 3t . 



4 SPHERICAL COLLAPSE 

Using the notation int roduced by IWang & SteinhardJ Jl998h and 
Bartelmann et al.l d2006t) . all relevant quantities are normalised by 
their values at turn-around, i.e. the scale factor a and the radius r of 
the density perturbation become 



a 

a„ 



y : 



(16) 



Generalising the Friedmann equations for incorporating the modi- 
fication of gravity on large scales outlined above yields (compare 
eqns.[4]and[3} 



and 



y 



y 
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ffl 



(17) 
(18) 



with to = Q.,„(a„) and A = I - to, assuring flatness. In this equa- 
tion, a perturbation with overdensity £ (measured at turn-around) 
will collapse at a specified time, f can be computed by solving the 
differential equation defining y(r) subject to the initial condition 
y(r = 0) = and the boundary condition y = 1 at turn-around, 
which is characterised by y = 0. At early times, g(x) =s x and the 
collapse equations reduce to the SCDM case, 



2y 2 ' 



(19) 



with the dot denoting derivatives with respect to the dimensionless 
time parameter r, 



t = H(a a )t. 



(20) 



The solution for x(t) can only be found numerically in the general 
casfl but simplifies to the SCDM result at early times x <k 1 or for 
A = 0: 

2 XT- 

T=~-i=- (2D 

i -yto 

The integration of y involves multiplying both sides with yy and 
noticing that 

1 _d_ 

2dr' 



>ys(ff- if™ (7 



(22) 



, y )) \ y 1 2 y 2 

which allows the reduction of the second-order differential equation 
for y to a first-order equation, 



dr 



, 2 



Subsequent integration yields 



+ c, 



(23) 



(24) 



where the integration constant is determined by requiring y = at 
turn around y = 1 : 

c = -g(a>0, (25) 
yielding the final result for the differential equation defining y(r): 



y 



y 



to£_ 



(26) 



For small radii y at early times t, the differential equation can be 
approximated by 



dr : 



due to g(f) a f at early times, and can be integrated to yield 
1 2 



1 + — — —y 



(27) 



(28) 



Vk3' \~ ' 10 

Comparing eqn. d28t with eqn. d2U at early times and eliminating 
t yields a solution for y(x), which can be used for defining the time 
evolution A(x) of the density f, restricted to early times: 



AW = (- 



y 



1 + 7 — y y- 

5 to£ 



(29) 



The linear density contrast S c is related to the overdensity A - 1 by 
computing the density resulting from linear growth, 



5 C = D + (x c ) lim 



M(x) 



- 1 



(30) 



x-K)\ £»+(x) 

The overdensity 6 C is shown in Fig. [5] as a function of redshift z, 
along with the constant value 6 C = 1.686 for spherical collapse 
in a SCDM cosmology. 6 c (z) falls significantly below the SCDM 
value at low redshifts with increasing a, reaching values as low 
6 C = 1.46 for modified gravity. Interestingly, the range of r5 r -values 
correspo nds to those obtained in early dark energy mode ls, as pro- 
posed bv lBartelmann et al.l d2006h and lSadeh etall ( 120071) . albeit at 
lower redshifts. Consequently, the observational signature of low- 
ered S c in modified gravity would be an increased abundance of 
clusters at the relevant redshifts, which will be probed especially 
by future Sunyaev-Zel'dovich surveys. 

1 An exact solution for the ACDM case (a = 0) can be derived in terms of 
the hypergeometric function, r(x) = ^S^l^lik 1 2' §' ~^j x3 ^ 
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redshift z 

Figure 5. Overdensity for spherical collapse 6 c (z), for a = (corresponding 
to ACDM, dash-dotted), a = 1 (DGP-gravity, solid) and intermediate values 
in steps of Act = 0.2 (dashed). Additionally, the constant value of 5 C = 1.686 
is plotted for spherical collapse in the SCDM cosmology. 
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b 


c 


a 


= 


1.6823 


0.0129 


1.5648 


a 


= 0.2 


1.6820 


0.0314 


1.0912 


a 


= 0.4 


1.6806 


0.0582 


0.9994 


a 


= 0.6 


1.6774 


0.0942 


0.9369 


a 


= 0.8 


1.6710 


0.1397 


0.8722 


a 


= 1 


1.6585 


0.1945 


0.8020 



Table 1. Fitting values for the empirical model S c (z) = a - foexp(-cz) for 
the range of the modified gravity parameter a considered here, which repro- 
duces the numerical values for <$ c (z) to sub-percent accuracy in the redshift 
range < z. < 5. 



The dependence of 6 C on redshift can be well fitted with the 
phenomenological function S c (z) = a - b exp(-cz). Table [T] sum- 
marises the values of the parameter a, b and c, which result from 
a fit to numerically obtained values for 6 c (z) in the redshift range 
< z < 5, for the range of modified gravity models (characterised 
by a) considered here. 

In the numerical solution, the sequence of steps is as fol- 
lows: From the specified collapse scale factor a c , the corresponding 
cosmic time t c is computed using t = fda/(aH) with the Hub- 
ble function H(a) obtained numerically from the Friedmann equa- 
tion eqn. (0. Due to the time-reversal symmetry of eqn. l |18t , turn 
around takes place at t a = t c /2, from which a a can be retrieved. 
The scale factor at turn around a a is used to define oj and A using 
eqn. H5\ . f is adjusted numerically by evolving eqn. d24t to the 
time t' at which y(r') = 0, and applying a numerical root finding 
scheme on t' - t(x =1), where the dimensionless time t(x = 1) at 
turn around x = 1 results from eqn. dl7t for given density param- 
eters oj and A. From these results, A(x) and S c (z) follow directly. 
In the numerical treatment, we use a very good approximation to 
the implicitly defined Hubble function H(a), which is described in 
detail in Appendix Icl 

It is worth noting that 6 C depends on the matter density Q.„, 
from the determination of ( via the turn around point in the evolu- 
tion of the ^-differential equation, as well as from the growth func- 
tion £> + . In the remainder of the paper, we neglect the influence of 



Q,„ on 6 C , as it is of minor influence regarding e.g. cluster number 
counts compared to first order dependence of the mass function on 

n,„. 



5 CLUSTER NUMBER COUNTS AND MERGER RATES 

Fluctuations in the cosmic density field on the scale k are described 
by the power spectrum P(k), which results from ensemble averag- 
ing of Fourier modes S(k), (6(k)6(k')) = (2nfS D (k + k')P(k). The 
shape of the power spectrum P(k) oc k" s T 2 (k) is wel l appro ximated 
by the transfer function suggested by Bardeen et al. ( 1986) 



T(q) 



ln(l +2.34q) 
2.34q 



(l + 3.89? + (\6.lq) 2 + (5.46q) 3 + (6.71^) 4 ) 5 (31) 



where the the wave vector k = qT is rescaled with the with shape 
parameter T ^ £l m h. The fluctuation amplitude is normalised to the 
value <r 8 on the scale R = 8 Mpc/h with a spherical top-hat filter. 

The number density of clusters forming by spherical collapse 
in a fluctuating density field is given by 

[2 d In cr(M) / A 2 (z,M)\dM 



dPress & Schechtedl 19741 ; iBond et al.lll99lh with background den- 
sity po = fl m p cr j t , where 

SAz) 



A(z, M) : 



D + (z)cr(M) 



(33) 



abbreviates the quotient between the linear overdensity S c (z) 
needed for collapse at z, the growth function D + (z) and the vari- 
ance <r(M) of the density field on the mass scale M, which is related 
to the spatial filtering scale R by M = 4w/3p cr i,f2 M A 3 . The Press- 
Schechter mass function is used rather than the extension by Sheth 
and Tormen because the latter describes ellipsoidal collapse, where 
the modified gravity collapse equations might bot be applicable. 

The logarithmic derivative of the fluctuation amplitude with 
logarithmic mass is reformulated as a derivative with respect to spa- 
tial scale R: 



diner 1 diner 



dlnM 3dlnfl 



R d<r 
3^^ 



R 2 dcr 2 

6cr 2 dR ' 



(34) 



where the fluctuation amplitude <r 2 (R) and its derivative d<r 2 /dr? 
are given by: 



<y\R) 

da 2 
~dR 



2n 2 J 



dktPWl—j^kR) 



dkeP(k)\ — \ h{kR)j 2 {kR). 



(35) 
(36) 



ji (x) denotes the spherical Bess el function of the first kind of order 
t iAbramowitz & Stegunll 19721) . 

Multiplication of n(M, z) with the evolution dV/dz of the co- 
moving volume with redshift yields the cluster density per unit red- 
shift, 



q(M, z) = n(M, z) 



dV 
d? 



(37) 



The quantity Q(M, z), which is defined as the ratio of q(M, z) for a 
modified gravity model characterised with a relative to ACDM, 



Q(M,z) 



q a (M,z) 

q a =o(M,z) 



(38) 
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Figure 6. Ratio Q(M, z) of the mass functions q(M, z) relative to ACDM 
(a = 0), at redshifts z. = 0.3 (thick lines) and z = 0.7 (thin lines), for values 
of a from a = 0.2 in steps of Ao- = 0.2 (dashed) up to ff = 1 (modified 
gravity, solid). 



Figure 7. Ratio R(M, z) of the merger densities r(m, z) relative to ACDM 
(a = 0), at redshifts z. = 0.3 (thick lines) and z, = 0.7 (thin lines), for values 
of a from a = 0.2 in steps of Ao- = 0.2 (dashed) up to a = 1 (modified 
gravity, solid). The mass ratio has been set to Mi/Mi = 2. 



is plotted in Fig. [6] for redshifts z = 0.3 and z = 0.7, where the 
PLANCK and SPT surveys are expected to yield most of their de- 
tections. 

Modified gravity significantly enhances the cluster number 
counts which is mainly caused by the lower values of the linear 
overdensity for spherical collapse 6 C at low redshifts. At a redshift 
of z = 0.3, corresponding to the peak in the redshift distribution of 
PLANCK's detections, on can expect to detect almost 50% more 
clusters at masses of l0 u M e /h. This effect is more pronounced at 
increased redshift: SPT is yield most of the detections at a redshift 
of z = 0.7, where number counts of 10 14 M o //z-mass objects would 
be almost twice as large in modified gravity models compared to 
ACDM. At higher masses of l0 [5 M o /h, the increase in cluster 
number density amounts to factors of three (at z = 0.3) up to six (at 
z=0.7), but the small spatial density of very massive objects at high 
redshift will make this difficult to observe. In contrast, the differ- 
ence between modified gravity cosmologies and ACDM vanishes 
at small masses, because low-mass objects form early the matter- 
dominated phase, where 6 C has not deviated significantly from its 
canonical value of 1.686. The enhancement in cluster number den- 
sity depends exponentially on the modified gravity parameter a. 

The density m(M, z) of merging events as a function of red- 
shift and mass ratio M 2 /M[ can be determi ned by extending 
the Press-Schechter formalism outlined above (Bond etal.lll99ll ; 
lLacev&Colej|l993l . ll994l) for an expression for the merger proba- 
bility as a function of mass, mass difference and cosmic time, 

d 2 p(M,z) 
dlnAM dt' 



(39) 



For the purpose of this work, the merger density is multiplied with 
the volume element dV/dz and the evolution dt/dz of cosmic time 
with redshift in order to yield the merger density per unit redshift, 



r(M,z) = \n(M,z) ■ 



dV\ ( d 2 p(M,z) dt 



dz \ d In AM dt dz 



(40) 



R(M, z) is defined as the ratio of r(M, z) between a modified gravity 
model and ACDM, 



R(M,z) 



r a (M,z) 
r a=0 (M,z)' 



(41) 



The quantity R(M,z) is plotted in Fig. [7] for redshifts z = 0.3 and 
z = 0.7 for a fixed merger mass ratio of M 2 /Mi = 2. Merger 
densities exhibit a similar behaviour at high masses of 10 l5 M o /h 
compared to the mass function, where the modified gravity model 
predics numbers roughly factors of 1.6 (PLANCK) and two (SPT) 
times higher at the relevant redshifts. Interestingly, the merger den- 
sity is decreased in modified gravity models relative to ACDM 
at small masses; this decrease is due to the unique weighting of 
the merger density with the evolution of cosmic time with redshift 
dt/dz, the evolution of comoving volume dV/dz and the Hubble 
function H(a), resulting from the conversion of a time derivative of 
S c /D + in eqn. l |39t to a derivative in scale factor. 



6 CONSTRAINTS FROM FUTURE SZ SURVEYS 

In this section, constraints on cosmological parameters includ- 
ing the modified gravity parameterisation a from future high- 
yield Sunyaev-Zel'dovich c luster surveys are es t imate d in a Fisher- 
matrix analysis. Following Majumdar & Mohr (2004), a measure- 
ment of the cluster number density per unit redshift interval, ob- 
served in the solid angle AQ., 



N(Z)1 



dn AQ. dV f° , w . 

— = — dM n(M, z), 

dz An dz J MnM (z) 



(42) 



is considered with a simple flux threshold S ra j„, expressed as a min- 
imal mass M min dependen t on redshift. Using the formula given by 
Maiumdar & Mohi] d2004l) . the flux distortion A5 (v) caused by the 
SZ-effect (in mjansky) at frequency v is related to the angular di- 
ameter distance d A (z) (in Mpc) and the cluster mass M (in M Q ) by 

AS(v) = -^-/icm gsz(v) A sz HHz) M 1 ^ (1 + z), (43) 

where /i CM = 0-12 a nd log A = -16.71 jDiego et all l200ll; 
iFinoguenov et alj|200ll) . For the slope of the temperat ure-mass re- 
lation , T(M) oc M?, the value p = 0.65 is chosen dMohr et all 
1999). The frequency dependence o f the Sunyaev-Zel'dovich effect 
dSunvaev & Zeldovichll 19721 . 1980) is given by 



x 4 exp(x) / exp(jc) + 1 



(exp(x) - l) 2 \ exp(x) - 1 



(44) 
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PLANCK 


SPT 


frequency 


v = 353 GHz 


v = 150 GHz 


flux 


S min = 100 mjy 


S mm = 10 mjy 


solid angle 


AH = 4tt 


An = 4jt/10 



Table 2. Frequency v, flux threshold 5 ra i n and the solid angle An of the 
Sunyaev-Zel'dovich surveys to be carried out by PLANCK and by SPT, 
respectively. 



with the dimensionless frequency x = hv/(k B T C MB)- In addition to 
the flux criterion, a lower mass threshold of l0 l4 M o /h is imposed, 
whichever yields the larger limiting mass. 

From the cluster number counts jV; in a redshift bin centered at 
Zj and the derivative dNj/dp^ with res pect to the parame ter p^ the 
Fisher matrix Fzr can be constructed (Hold er et al]200ll) : 



z 



dNj 1 dNj 
dp M Nt dp v ' 



(45) 



assuming a Gaussian likelihood and statistically independent num- 
ber counts Nj in each redshift bin.The summation extends over the 
redshift range from Zmm = 0.01 to z max = 2.0 in steps of Az = 0.02. 

Using the Cramer-Rao inequality, the accuracy on a single 
parameter Ap^ is equal to ^j(F' l ) m . x 2 for a pair of parameters 
(Pm>Py) can b e computed from the inverse (F -1 )^ of the Fisher 
matrix, 



(46) 



where Ap^ = p M — p^ CDM . For parameter pairs the isoprobability 
contour at 



■ In erfc 



(47) 



encloses a fraction of erf(n/ V2) for of the admissible parameter 
space, corresponding to a confidence level of ncr. erfc(x) = 1 — 
erf(x) is the complementary error function dAbramowitz & Stegun 
1 19721) . 

For the purpose of this paper, the PLANCK and SPT 
Sunyaev-Zel'dovich surveys are considered, which are expected 
to yield cluster c atalogues compris i ng 10 4 and 3 x 10 4 entries , 
respectively fc.f. iMelin et alj l2005t iGeisbusch & Hobsonl 120061 : 
ISchafer & Bartelmannl2007l) . The principal characteristics of these 
surveys are summarised in table [2] As cosmological parameters, 
Cl m (assuming spatial flatness), cr 8 and n s are considered, alongside 
the modified gravity parameter a and the slope of the temperature- 
mass relation /3. In general, the increased spatial number density of 
clusters caused by deviations in the gravitational law has to com- 
pete against a slower evolution of the comoving volume with red- 
shift dV/dz, since the two terms in eqn. J42b are affected in opposite 
ways. 

The results from the Fisher matrix analysis are compiled in 
Fig. [D The SZ-observations alone yield accuracies of a few per- 
cent on the parameters Q.„, and cr 8 , and marginally worse con- 
straints on n s . PLANCK and SPT would be able to exclude values 
of a > 0.4, which compare fa vourably to contemporary c onstraints 
from SNIa observations (e.g. iFairbairn & Goobaal 20061 who find 
the entire range of a consistent with observations). SPT's smaller 
power in constraining a is due to the fact that its SZ-sample con- 
tains more high-redshift clusters, which have formed in an epoch 
close to matter domination where the influence of modified grav- 



ity is smaller. An important consistency check is that the accuracy 
of measuring a derived from the SZ-surveys compares well to the 
allowed range of values fo r the dark energy equation of state w(a) 
dMaiumdar & Mohj|2004h of Aw = 0.1 (c.f. Fig.|ATJ where the 
w(a) for a dark energy model with the identical expansion history 
as the modified gravity model is derived). 

The a-crg degeneracy is readily explained by the fact that clus- 
ters form more easily in modified gravity, which has to be compen- 
sated by lowering the value of cr 8 . Tilting the CMB spectrum and 
using smaller values for n s yields more massive clusters if the nor- 
malisation tr 8 is kept fixed. Whether this is observable is a matter 
of the redshift distribution of the experiment: At low redshifts of 
z — 0.3, where most of the PLANCK-detections will be situated, 
the opposite effects of increasing a on the cluster density and on 
the volume element almost cancel each other, whereas at high red- 
shifts of z — 0.7, corresponding to the redshift at which SPT will 
find most of the objects, the decrease in dV/dz caused by modified 
gravity dominates. 

In the a-Q.„, degeneracy both parameters impact on the cluster 
density and the evolution of the comoving volume element: At low 
redshifts, the increased number densities of clusters due to modified 
gravity is balanced by low choices of Q m , which affect both the 
cluster number density as well as the volume element. At higher 
redshifts, however, the volume element becomes the dominating 
factor as it is decreased by modified gravity, only balanced by larger 
comoving densities due to higher choices of Sl m . 

A crucial point is the understanding of flux threshold due to 
the strong dependence of the observed number of clusters on the 
particular value of the limiting minimal mass as a function of red- 
shift. The key parameter is the slope of the temperature-mass re- 
lation p, which has to be known to the level of a few percent, 
and which would be cont rolled by self-calibration techniques (e.g. 
iMaiumdar & Mohj|2004l) in a realistic application. 

Additionally, Fig. [9] shows three-dimensional degeneracy el- 
lipsoids between a pair of cosmological parameters and the modi- 
fied gravity parameter a. It emphasises the fact that a strong prior 
on the cosmological parameters relevant for cluster formation can 
significantly improve the constraint on modified gravity. Especially 
SPT should be able to yield a constraint on a with independent 
measurments of £l m , n s and especially erg. 

Strong prior constraints on fl m , cr 8 and n s are given by the 
observation of e.g. CMB temperature and polarisation anisotropies: 
For independent observational channels, the combined likelihood 
of the parameters is the product of the individual likelihoods, and 
hence 



C-SZ , pCMB 

1 flV flV ' 



(48) 



where F™ B is derived for PLANCK's observation of tempera- 
ture anisotropy spectrum C T {£), the polarisation anisotropy spec- 
trum Ce(£) and the cross-spectrum Cc(C)- The construction of the 
Fisher matrix F™ B and details such as the assumed noised levels 
and beam dimensions are given in Appendix [D] 

Expected errors on individual cosmological parameters with 
and without CMB priors are given in table [3] Apart from percent 
accuracy on the cosmological parameters Q m , erg and n s , the min- 
imum variance bound on the modification of gravity derived from 
the Fisher-matrices amounts to 0.46 from PLANCK and to 0.29 
from SPT, i.e. from a cluster survey alone one should be able to 
exclude modified gravity models (corresponding to a = 1) with a 
statistical significance of 3<r. The improved measurement is able to 
constrain the modified gravity parameter a significantly better to 
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0.28 




0.8 1 1.2 




0.645 0.65 0.655 




Figure 8. Forecasts for parameter degeneracies for the PLANCK (light gray) and SPT (dark gray) Sunyaev-Zel'dovich surveys in a Fisher-matrix analysis. 
The fiducial ACDM cosmology is marked by a circle: fi m = 0.25 (assuming spatial flatness), erg = 0.8, n s = 1.0, the modified gravity parameterisation a = 
and the slope of the temperature-mass relation /J = 0.65. The ellipses correspond to \cr and 2cr confidence intervals. 




Figure 9. Simultaneous constraints on three cosmological parameters from the PLANCK (light gray) and SPT (dark gray) Sunyaev-Zel'dovich surveys, for 
the parameters Sl m , erg and a (left panel) and erg, n s and a (right panel). The circles mark the parameter choices fi„, = 0.25, erg = 0.8, n s = 1 and a = for 
the fiducial ACDM cosmology. The nested ellipsoid surfaces correspond to \cr and 2cr confidence intervals. 
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PLANCK 


SPT 




PLANCK + CMB 


SPT + CMB 


An,,, = 0.027 


Afi,„ 


= 0.008 


AD. m = 0.002 


AQ m = 0.001 


Ao- 8 = 0.082 


Ao- 8 


= 0.021 


Ao- 8 = 0.002 


Ao-8 = 0.002 


An s = 0.023 


An s = 


= 0.057 


An s = 0.002 


An s = 0.002 


Aa = 0.469 


Aa = 


0.292 


Aa = 0.063 


Aa = 0.035 



Table 3. Fisher-matrix forecasts for parameter accuracies (lo - confidence 
intervals) from the PLANCK and SPT Sunyaev-Zel'dovich surveys (first 
and second column), and the corresponding accuracies after adding prior 
constraints on Cl m , o"8 and n s from PLANCK's CMB temperature and polar- 
isation spectra (third and fourth column). ACDM is the fiducial cosmology, 
with parameters Sl m = 0.25, o"s = 0.8, n s = 1 and a = 0. 



values below Aff = 0.063 and Aa = 0.035, for PLANCK and SPT, 
respectively. 

These resul ts are quite comparable to the findings of 
iTang et alj l l2006h . who combined supernovae and cluster counts 
for constraining DGP-gravity, which is achievable with a few hun- 
dred clusters from an SZ-observation, given a prior tight constraint 
on cr 8 . In their analysis, the influence of modified gravity on the 
Hubble expansion and on the growth equation was considered, and 
the density t hreshold by spherical collapse was approximated to 
be constant. ITang et alT d2006h point out that a pure geometrical 
measurement can never distinguish between modified gravity and a 
suitably tuned dark energy model, and that additional information 
such as a measurement of the growth function is needed. 



7 SUMMARY 

In this paper, the influence of modified gravity on Friedmann dy- 
namics, the homogeneous growth of the cosmic density field and 
the formation of galaxy clusters is considered. The alteration of the 
gravitational law is motivated by DGP-gravity, and introduced as an 
additional term in the Friedmann equations, which are consistently 
applied to the background dynamics as well as for the formation of 
individual clusters. A parameter a allows interpolation between the 
two limiting cases: gravitation with a cosmological constant A and 
DGP-gravity. 

• The Friedmann dynamics of the universe is remarkably little 
affected by the modification of gravity, as the Hubble function H(a) 
for our range of models differs by at most 10% at a scale factor of 
a = 0.5. This translates to comparable variations on the evolution 
dV/dz of volume with redshift and the evolution dt/dz of cosmic 
time with redshift. 

• The growth function shows significant suppression of growth 
at early times for modified gravity, very similar to dark energy mod- 
els with evolving equation of state, where the first derivative of the 
equation of state parameter w{a) is positive at zero redshift. The 
largest spread between the models considered here occurs around 
a = 0.3. The changed behaviour of the growth function is both 
caused by the decrease in matter density Cl„,(a) as well as by the 
scaling of the Hubble-function with scale factor. 

• The spherical collapse equations can be solved numerically 
for the modified gravity models motivated by DGP, resulting in the 
linear overdensity for spherical collapse S c as a function of collapse 
redshift z. One notices a significant decrease in S c by almost 20% 
compared to SCDM. The range of values for S c is quite comparable 
to those obtained in collapse with early dark energy. An empirical 



fitting formula of the shape S c (z) = a - bexp(-cz) was found to 
yield a good approximation. 

• Cluster number counts are significantly enhanced by modified 
gravity, due to the lower values of the linear overdensity for spher- 
ical collapse 8 C at low redshifts. At a redshift of z = 0.3, where the 
majority of PLANCK SZ-clusters will be situated, as much as 50% 
more clusters with masses of l0 14 M Q /h have formed, whereas at a 
redshift of z = 0.7, corresponding to the peak in the cluster redshift 
distribution forecast for SPT, number counts might be 70% larger. 
At higher masses of 10 I5 M o //i, the increase in cluster number den- 
sity amounts to factors of three (at z = 0.3) up to six (at z=0.7), 
but the sparsity of very massive objects at high redshift will make 
this difficult to observe. At low masses, on the contrary, there is 
hardly a difference between DGP-like modified gravity models and 
ACDM, which is due to the fact that these objects form early well 
in the matter-dominated phase. The enhancement in cluster number 
density depends exponentially on the modified gravity parameter a. 

• Merger densities exhibit a similar behaviour at high masses 
of l0 i5 M o /h, where the model predics an increase by factors of 
roughly up to two at the redshifts considered. At low masses, how- 
ever, the merger density falls below the prediction of ACDM, which 
is carried by the unique weighting of the merger density with the 
evolution of cosmic time with redshift dt/dz, the evolution of co- 
moving volume dV/dz and the Hubble function H(a). 

• The impact of modified gravity on cluster formation is most 
notable at small redshifts below unity, making it a target for 
Sunyaev-Zerdovich cluster surveys. In this paper the cluster yield 
of the SZ-surveys to be carried out by PLANCK and SPT was 
modelled with a simple flux threshold as the detection criterion. 
The Fisher-matrix analysis with ACDM as the fiducial cosmologi- 
cal model found that on themselves, the SZ-surveys alone are able 
to constrain the modified gravity parameter a with an accuracy of 
Aa =! 0.47, albeit with strong degeneracies. Adding in CMB con- 
straints on Q.,„, erg and ns significantly lowered the uncertainty to 
Aa = 0.04, making the combined SZ- and CMB-observations a 
powerful observational probe. 

As last points is it worth noting that the above described 
model for phenomenological DGP-gravity has the same number of 
parameters as ACDM and would not be disfavoured on grounds 
on Bayesian model selection. The feature that cluster formation 
takes place at higher redshifts compared to ACDM might have 
imp ortant implications on t he arclet problem in gravit ational lens- 
ing (Meneg hetti et aljl 20051) and on the CBI-anomaly (Pa din et"al] 
12001 ; Bond et al. 20051) as it naturally causes a higher level of 
Sunyaev-Zel'dovich background due to high-redshift clusters with 
a low contemporary value of <x 8 . 
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APPENDIX A: EQUATION OF STATE 

The relation between the modified gravity model and the equation 
of state w(a) of a dark energy model with identical expansion his- 
tory is derived in this section. A dark energy fluid with equation of 
state w(a) gives rise to an additional term AH 2 (a) in the Friedmann 
equation, given by 



AH (a) = exp 



1 1 + w(a) 
da 



(Al) 



which can be solved for w(a) for a given Hubble function H(a). By 
equating AH 2 {a) with the modification in the Friedmann equation 
( 1 - Q.,„)H a (a) the equation of state vv^ (a) can be derived for a dark 
energy model with identical expansion history: 



o-dlntf 
1 + 

3 din a 



(A2) 



Another interesting quantity is the equation of state w avg (a), av- 
eraged over all relevant cosmological fluids. In this case, AH 2 
has to be equal to the complete Friedmann equation H 2 (a) = 
(1 - fi,„ )H a (a) + £l m /a 3 , which results in 



w„Ja) : 



1 + 



2a din H 
T dlna 



(A3) 



Fig. lAll compares the average equation of state w mg (a) and the ef- 
fective equation of state w cff (a) as a function of scale factor. The 
effective equation of state w c ff(a) decreases slowly from -0.5 to - 
0.8 for the DGP-type models, is constant w = -1 for ACDM and 
smoothly interpolates between these two limiting cases for choices 
of a = 0. . .1. On the contrary, the average equation of state w avg (a) 
assumes a value of zero at a = for all modified gravity mod- 
els, and decreases towards values in the interval -0.6. . . - 0.8, with 
ACDM showing the most abrupt and the DGP-emulating model 
showing the most gradual transition. 

Comparing the results in Fig. lAll to current measurements of 
the equa tion of state of dar k energy, e.g. using weak gravitational 
lensing |jarvis et ail 12006)) . th e cosmic m icrowave background 
dVerde et al.l2003l)7supemo vae jRiessl20 04). baryon acoustic oscil- 
lations ( Percival et al .l2007h or combinations of these observational 



channels jBean e t al. 200 1 ; Melchiorri et al.ll2003l) , which are usu- 
ally stated in terms of the parameterisation w(a) ^ wq + (1 - a)w a 
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Figure Al. Average equation-of-state parameter w avg (a) (equal to at early 
times) and effective equation-of-state parameter w e g(a) as a function of 
scale factor a, for values of a = (corresponding to ACDM, dash-dotted), 
a = 1 (DGP-gravity, solid) and intermediate values in steps of Aa = 0.2 
(dashed). 
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Figure Bl. Comparison of the damping term Q(a) = 3 + dln//(a)/dlna 
(equal to 3/2 at early times) and the source term S(a) = G(a)Sl m (a) (equal 
to 1 at early times) in the growth equation as a function of scale factor a, for 
values of a = (corresponding to ACDM, dash-dotted), a = 1 (DGP-like 
gravity, solid) and intermediate values in steps of Aa = 0.2 (dashed). 



(introduced by lLinder& Jenkins! 120031) shows that these experi- 
ments only marginally exclude DGP-gravity and are in agreement 
with the values of wq = -0.9. . . - 1.0 and w a = 0.3. . .0 resulting 
from more conservative choices of a. 



APPENDIX B: GROWTH FUNCTION 

The influence of modified gravity on the growth equation is dis- 
cussed in this section. Generally, The growth equation eqn. (13t 
has the shape 



d 2 Q(a) d 35(a) 

—D + {a) + ^—D + (a) = D + (o), 
aa z a aa la 1 



(Bl) 



where a source term 5 (a) and a dissipation term Q(a) can be iden- 
tified, 



Q(a) = 3 + 



din H(a) 
din a 



5(a) = G(a) Q m (a)- 



(B2) 



A comparison between the source and damping terms in the growth 
equation is provided by Fig. IB 11 Increasing the value of a makes 
the transition of the damping term Q(a) from Q(a = 0) = | in the 
matter dominated epoch to the modified gravity dominated epoch 
more gradual, due to the influence of modified gravity on the scal- 
ing of the Hubble function. 

The source term 5(a) = G(a) Q.„,(a) with contributions both 
from the evolution fl,„ of the matter density, which is proportional 
to H~ 2 (a) as well a s from modified gravity as a consequence of 
Birkhoff 's theorem fcue et al.|[2004al lbh. smoothly decreases from 
5 (a) = 1 at a = to 5 (a) = Q.,,, today, and the choice of larger val- 
ues of a make the source term drop at earlier times. In summary, 
the suppression of structure formation is caused simultaneously by 
the decrease in the source term 5 (a) and the increase in the damp- 
ing Q(a), while larger choices of a are able to suppress structure 
formation more effectively at early times. 



APPENDIX C: APPROXIMATE SOLUTION FOR THE 
HUBBLE-FUNCTION 

As discussed in Sect. [3] our family of modified Friedmann equa- 
tions, 

H\a) = (1 - fi,„)tt» + (CI) 
a J 

can only be solved numerically for H(a). In order to have a much 
faster analytic approximation for the spherical collapse equations, 
we use two second order expansions, yielding quadratic equations, 
which are analytically solvable. At late times a = 1, we expand 
H a (a) at H = 1 to second order, 

a(a - 1) , 
H" = 1 +a(H- 1) + — - — -{H-lf, (C2) 

substitute into eqn. JC 1 fr . and solve the resulting quadratic equation 
for the approximation H + (a), which is valid at late times. At early 
times, where the Hubble function diverges, we expand Ho(a) = 
a?H(a), which assumes the value y/£l„, at a = 0: 

H a - nj + aQ.f(H - VO^) + a(a ~ l) Slf{H - VOJ 2 . (C3) 

This expansion can be converted to a second order expansion for 
H(a), substituted into eqn. dCU and the resulting quadratic equation 
solved for the early time approximation H-(a). The relation 



H{a) = y_ (a) H_(a) + y+(a) H + (a) 



(C4) 



interpolates between the two branches by adjusting the weight a 
accordingly, and avoids any sharp transition. y ± (a) is a nonlinear 
interpolation, 



1 

r±(«) = J 



1 ± — arctan 



Pi 2 ) 



(C5) 



with the transition scale factor fi = 0.3 and the transition width 
A = 0.05. Eqn. dC5t is normalised in the sense y + (a) + y_(a) = 1. 
Figure |CT| compares the approximation of the Hubble function to 
the direct numerical solution, and gives the relative deviation be- 
tween the two: The approximation reproduces the Hubble function 
to an accuracy better than 0.3%, which is accurate enough for the 
purpose of spherical collapse computations. It should be noted, that 
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Figure CI. The scaled Hubble function H(a) from the numerical approxi- 
mation, as a function scale factor a, for values of a = (corresponding to 
ACDM, dash-dotted), a = 1 (DGP-gravity, solid) and intermediate values 
in steps of Aa = 0.2 (dashed). At the bottom, the relative deviation in per- 
cent between the direct numerical solution for H(a) and the approximation 
is plotted. 



and Ct(£) = Cc(£)- For PLANCK's noise levels the values w T l = 
(0.02p_K) 2 and Wp = (0.03 p.K) 2 have been used, and the beam was 
assumed to be Gaussian, B~\£) = exp(A8 2 £(£+ 1)), with a FWHM- 
width of Ad = T.l, corresponding to the v = 143 GHz channel 
closest to the CMB -maximum. The observation uses a fraction of 
/sky = 0.8 of the entire sky, corresponding to a galactic cut exclud- 
ing \b\ < 11.5 °. CMB spectra an d their derivatives were generated 
with CAMB0 dLewis et alfcOOOl) . 

This paper has been typeset from a TgX/ KTgX file prepared by the 
author. 



the cases a = (ACDM) and a = 1 (DGP-gravity) are reproduced 
exactly. 



APPENDIX D: ADDITIONAL CONSTRAINTS FROM 
CMB TEMPERATURE AND POLARISATION SPECTRA 

The constraints on modified gravity from cluster number counts 
can be significantly improved by adding priors on Q.,,,, erg and n s 
from CMB temperature and polarisation spectra. The Fisher matrix 
of the CMB temperature and polarisation power spectra is given by 
JZaldarriaga eTatll 19971 : iBalbi et alj|2003l) : 



= X X -^Cov- 1 [C X (£)C Y ({)] 
e=2 x,y " x '< " Xv 



(Dl) 



where the indices (X, Y) assume values of T for the temperature 
spectrum Cr(i), E for the £-mode polarisation spectrum Ce{() and 
C for the temperature-polarisation cross-spectrum Cc(£)- 

The entries of the covariance matrix Cov [C X (£)C Y (£)] are 
given by: 

2 1 



Cov [CA£)C T (0] 
Cov [C E (£)C E (£)] 
Cov [C C (£)C C (£)] 
Cov [C T (l)C E (m 
Cov [C E (£)C c (m 
Cov [C T (£)Cc(£)] 



C 2 T (£) 



C\{£) 



2£+lf sky 
2 1 

2£+ l]ky 

2 1 

2£+l ]ky 
-1 -C 2 C (£), 

2 1 



2^+1 / Sk y 

2 1 

2£+ lJTy 



(D2) 
(D3) 

(C 2 C (£) + C T (£)C E (£)) , (D4) 
(D5) 

C E (£)C C (£), and (D6) 



C T (£)C C (£), 



with the observed spectra 

C T (£) = C T {€) + Wj l B~\€), 
C E (£) = C e (£) + w-;b- 2 (£), 



(D7) 



(D8) 
(D9) 



http://camb.info/ 



